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PACS 32.10. Fn - Fine and hyperfine structure 
PACS 32.70. Jz - Line shapes, widths, and shifts 

Abstract - A one-dimensional nano-metric-thin cell (NC) filled with potassium metal has been 
built and used to study optical atomic transitions in external magnetic fields. These studies 
benefit from the remarkable features of the NC allowing one to use A/2- and A-methods for 
effective investigations of individual transitions of the K D\ line. The methods are based on 
strong narrowing of the absorption spectrum of the atomic column of thickness L equal to A/2 
and to A (with A = 770nm being the resonant laser radiation wavelength). In particular, for 
a TT-polarized radiation excitation the A-method allows us to resolve eight atomic transitions 
(in two groups of four atomic transitions) and to reveal two remarkable transitions that we call 
Guiding Transitions (GT). The probabilities of all other transitions inside the group (as well as the 
frequency slope versus magnetic field) tend to the probability and to the slope of GT. Note that 
for circular polarization there is one group of four transitions and GT do not exist. Among eight 
transitions there are also two transitions (forbidden for B — 0) with the probabilities undergoing 
strong modification under the influence of magnetic fields. Practically the complete hyperfine 
Paschen-Back regime is observed at relatively low (~ 1 kG) magnetic fields. Note that for K D 2 
line GT are absent. Theoretical models describe the experiment very well. 
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Introduction. — Atomic spectroscopy with NC filled 
O with Rb or Cs atomic vapor, with a thickness of the vapor 
column L which is of the order of optical radiation wave- 
CN length has been found to be very efficient to study optical 
atomic transitions in external magnetic fields [T]. There 
are two interconnected effects: splitting of the atomic en- 
. 'ergy levels to Zeeman sublevels and shifting of frequen- 
. ^ (deviating from the linear dependence observed in 
quite moderate magnetic fields), and significant change 
in atomic transitions probabilities as a function of the B 
Afield PHI]. These studies benefit from the following fea¬ 
tures of NC: i) sub-Doppler spectral resolution for atomic 
vapor thickness L = A/2 and L = A (A being the resonant 
wavelength of Rb Di ^2 or Cs £> 1^2 lines) needed to resolve 
a huge number of Zeeman transition components in trans¬ 
mission or fluorescence spectra; ii) possibility to apply a 
strong magnetic field using permanent magnets in spite of 
a strong inhomogeneity of the B field (in our case it can 
reach 150G/mm). Note that the variation of the B field 
inside the atomic vapor is negligible as the vapor column 
thickness is small. K vapor in magnetic helds was stud¬ 


ied in [5] using a few centimeter-long cell and saturation 
absorption technique. However, due to the presence of 
strong cross-over resonances in the absorption spectrum, 
the technique is useful only for B < 10 G. 

Here we report the first studies of K vapor confined 
in NC and under the influence of relatively low magnetic 
fields {B < 2kG). The NC with L = A/2 and A are 
used for (T+ and tt laser polarizations. Also two theoret¬ 
ical considerations are taken into account. The splitting 
of atomic levels in weak magnetic helds is described by 
the total angular momentum F = J + I of the atom and 
its projection mp, where J = L + S is the total angu¬ 
lar momentum of electrons and I is the nuclear spin. In 
the hyperhne Paschen-Back (HPB) regime, J and I be¬ 
come decoupled and the splitting of the atomic levels is 
described by the projections mj and mj. For alkali met¬ 
als the hyperhne Paschen-Back regime takes place at helds 
B ^ Bq = Ahfs/p-B, where Ahfs is the ground-state hy¬ 
perhne coupling coefficient and is the Bohr magneton. 
For ^^^Cs, ®^Rb, and ®^Rb atoms Bq 1.7kG,2kG and 
0.7 kG, respectively isHinj- It is worth noting that 
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Fig. 1: The photograph of the K NC with a tapered gap of pjg 2: The hyperfine (hfs) energy levels diagram of the Di line 
thickness L ranging from 50 to 1500 nm. TC is a thermocouple of (the prime is used for the upper levels), 
to measure the side-arm (SA) temperature. 


has the smallest Bq « 0.17kG caused by the smallest hy¬ 
perfine splitting A Ri 462 MHz of the ground level 4 S'i /2 
[note that Ahfs = ^/(^ + 1/2)]- Thus, one expects to 
obtain complete J and I decoupling (HPB) of at rela¬ 
tively low magnetic fields B ^ 170 G. The manifestations 
of HPB regime, particularly, are as follows: i) strong re¬ 
duction of the number of atomic transitions (compared 
with that at low magnetic held) to a hxed number which 
is easy to determine from the diagram on the basis of mj 
and mi projections, ii) the hxed frequency slope inside the 
group of transitions, iii) the transition probabilities tend 
to the same value inside the group, iv) the energy of the 
ground and upper levels can be calculated from analytic 
expression &■ 

Experiment. — For our purpose a one-dimensional 
nano-metric-thin cell hlled with natural potassium (93.3% 
^®K, and 6.7% ^^K) has been built (for the hrst time) and 
used for the experiment. The design of the NC is similar 
to that of the extremely thin cell described earlier [TT] (the 
details of NC design can be found in [H]). The NC allows 
one to exploit a variable vapor column thickness L in the 
range of 50 — 1500 nm. It is demonstrated that for K vapor 
the key parameter determining the spectral width and the 
shape of the absorption line in the NC is the ratio A/A, 
with A = 770 nm being the wavelength of the laser radi¬ 
ation resonant with the atomic transition of Di line (as 
it was earlier demonstrated for the NC filled with Cs or 
Rb). In particular, it was shown that the spectral width 
of the resonant absorption reaches its minimum value at 
L = (2n -I- l)A/2 (n is an integer); this effect has been 
called the Dicke-type coherent narrowing [nmnn]. It is 
also demonstrated that for L = nX the spectral width of 
the resonant absorption reaches its maximum value, close 
to the Doppler width (> 0.9 GHz). 

A photograph of the K NC with a tapered gap is shown 
in Fig. [T] The windows of the NC were constructed with 
well-polished crystalline sapphire with the c axis perpen¬ 
dicular to the window surface to minimize birefringence. 
The regions with the thickness A = A/2 = 385 nm and 
A = A = 770 nm are marked by ovals in Fig.[l] A sapphire 
side-arm (SA) filled with metallic K is seen at the bottom 
of Fig. [T] The SA was heated to 150 — 160 °C (the temper¬ 
ature on the windows was by 20 degrees higher in order 



Fig. 3: Sketch of the experimental setup. ECDL, diode laser; 
FI, Faraday isolator; 1, A/4 plate; main K NC with thickness 
A = A/2 or A = A in the oven; PBS, polarizing beam split¬ 
ter; 2, auxiliary K NC; 3, 3' permanent magnets: in case of 
radiation use, the B-field is directed along the laser propaga¬ 
tion direction k (magnets 3); for tt polarization, the B-field is 
directed along the laser electric field E (magnets 3'); 4, pho¬ 
todetectors; 5, oscilloscope. 




(a) (b) 

Fig. 4: Absorption spectrum of vapor contained in the 
NCA = A/2 = 385 nm a) experiment, the laser power is 
2 /iW, the SA temperature is 150 °C, a small absorption of 
is also seen; b) theory [16] with the parameters: thermal 
velocity vth = 450 ms“^, Rabi frequency Q./2'k = O.OGyw (7iv ~ 
6 MHz). 

to prevent vapor condensation) providing the density of 
atoms A = (5 — 8) X 10^^ cm“^. The NC can operate up 
to temperatures of 500 °C. The hfs energy levels diagram 
of the Di line of is shown in Fig. [2] A sketch of the ex¬ 
perimental setup is shown in Fig.|3l The linearly polarized 
beam of an extended cavity diode laser ( 7 ^ < 1 MHz), res¬ 
onant with a Di line after passing through a Faraday 
isolator (FI), was focused onto a 0.5 mm diameter spot on 
the K NC orthogonally to the cell window. A PBS was 
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Fig. 5: Four VSOPs located at the positions of the transitions 
1, 2 —>■ 1', 2' are presented in the transmission spectrum of 
vapor contained in the NC with L = 770 nm; a small absorp¬ 
tion of is also seen; the SA temperature is 160 °C. The in¬ 
set shows the results of fitting by Lorentzian profiles; note that 
the ratio of the VSOPs amplitudes coincides with that deduced 
from the values given in Fig. [2] which means linear dependence 
on laser power. The lower curve - the theory [1^ for the pa¬ 
rameters: Vth = 450ms“^, Rabi frequency Q/2tv — 0.2yN- 


used to purify the initial linear polarization of the laser; a 
A/4 plate (1) was used to produce a circular polarization. 
In the experiments the thicknesses of the vapor column 
L = X and L = X/2 were used. The transmission signal 
was detected by a photodiode (4) and was recorded by 
a Tektronix TDS 2014B four-channel storage oscilloscope 
(5). To record the transmission spectra, the laser radia¬ 
tion was linearly scanned within up to a ~ 5 GHz spectral 
region covering the studied group of transitions. About 
30% of the laser power was branched to the reference unit 
with an auxiliary K NC (2). The absorption spectrum of 
the latter with thicknesses L = X/2 or L = X was used as 
frequency reference. 

The experimental and theoretical absorption spectra of 
K atomic vapor contained in the NC with L = X/2 = 
385 nm are shown in Fig. H] (the large width of the spec¬ 
trum at the base is caused by huge Doppler-broadening 
> 0.9 GHz at 170 °C at the NC windows). There is 
a strong spectral narrowing which allows one to sepa¬ 
rate four transitions between hyperfine sublevels shown 
in Fig. m The sharp (nearly Gaussian) absorption near 
the top makes it convenient to separate closely spaced in¬ 
dividual atomic transitions in an external magnetic field 
(we called it A/2-method [10] )■ We have also used the so 
called A-method [ISIIISI. In this case spectrally narrow 
velocity selective optical pumping (VSOP) resonances lo¬ 
cated exactly at the positions of atomic transitions ap¬ 
pear in the transmission spectrum of the NC with thick¬ 
ness L = X shown in Fig. [0| The VSOP parameters are 
shown to be immune against 10% thickness deviation from 
L = X, which makes the A-method feasible. Laser power 
is 15 ^W(for a lower power the VSOP line-width is less 


Fig. 6: Absorption spectrum of vapor contained in the NC 
with L = X/2 for B — 400,700 and 1440 G; excitation. 
The bottom curve is the reference showing the positions of ®®K 
transitions for B — 0. The inset shows the results of fitting 
by four ” pseudo-Vigt” functions (note that the amplitude of 
the transition labeled 4 is larger than that of 3). The absolute 
value of the peak absorption of transition 1 is 0.3%. 

than 30 MHz, but is more noisy). In a magnetic field, the 
VSOPs are split into many components. The amplitude 
and frequency positions of VSOPs depend on the H-field, 
which makes it convenient to study each individual tran¬ 
sition separately [1]. 

Magnetometry with A/2- and A-methods, cr+ 
laser excitation. — The assembly of an oven with NC 
inside was placed between the permanent ring magnets. 
The magnetic field B was directed along the laser radiation 
propagation direction k. An extremely small thickness 
of the NC is advantageous for the application of strong 
magnetic fields with the use of permanent magnets (hav¬ 
ing a 2mm diameter hole for the laser beam passage). 
Such magnets are unusable for centimeter-long cells be¬ 
cause of the strong inhomogeneity of the magnetic field, 
while in NC, the variation of the H-field inside the atomic 
vapor column is by several orders less than the applied 
B value. The H-field strength was measured by a cali¬ 
brated Hall gauge with an absolute imprecision less than 
30 G throughout the applied B-field range. The absorp¬ 
tion and transmission spectra of vapor contained in 
the NC with L = X/2 and L = X (for tT+ radiation) versus 
magnetic field are shown in Fig. IHjand Fig. [T] respectively 
(the bottom curves are the reference ones). Four atomic 
transitions (as predicted for the HPB regime by the dia¬ 
gram presented in Fig. [5|) are well observable. Note that 
the frequency separation between the two well resolved 
VSOP labeled 1 and 2 is about 90 MHz, which is more 
than 10 times less than the Doppler line-width. Also, the 
VSOPs have the same amplitudes which is an additional 
evidence of HPB regime. Let us briefly stress the bene¬ 
fits of A/2- and A-methods. The A/2-method requires less 
laser intensity and in case of small absorption the peak 
absorption A is proportional to aNL, where a is the ab- 
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Fig. 7: Transmission spectrum of vapor contained in the 
NC with L = X for B = 230, 400, 572 and 917 G and excita¬ 
tion. The bottom curve is the reference showing the positions 
of transitions for B — 0. 

sorption cross section proportional to d? (with d being 
the dipole moment). Thus, by direct comparison of Ai 
(peak amplitudes of the absorption of the f-th transition), 
it is straightforward to estimate the relative probabilities 
(i.e. line intensities). On the other hand, the A-method 
provides a fivefold better spectral resolution. Thus, the 
methods can be considered as complementary depending 
on particular requirements. Moreover, it is easy to switch 
from A/2 to A just by vertical translation of the NC. Note, 
that in the case where there is a big frequency separation 
between the transitions (as it is for the ®’’Rb isotope) also 
a 1 mm long cell can be used mm- 

Theoretical simulations for the B-field dependences of 
the atomic transition frequency shifts and relative tran¬ 
sition probabilities for 1,2 —?► l',2' transitions of the K 
Di line were based on the calculation of the eigenvalues 
and eigenvectors of the Hamilton matrix versus magnetic 
field for the full hyperfine structure manifold HHHH]. Al¬ 
though for the K Di line the Rabi-Breit formula could 
be also used, however the theoretical model is preferable 
since it is valid also for K D 2 line. Note, that in the case 
of HPB regime the energy of the ground -^ 81/2 and upper 
4 Pi /2 levels for the Di line is given by the following 
formula [^14]: 

E\jmjimi) = Ahfsmjmi -b fiB{gjmj + gimi)B. (1) 

The values for the fine structure {gj) and the nuclear 
{gi) Lande factors and the hyperfine constants Ahfs are 
given in m- The dependence of the frequency shifts on 
the magnetic field (relative to the position of the 2 ^ 2 ' 
transition at R = 0) is shown in Fig. 9, where the red 
lines are obtained by exact numerical calculation, while 
the black lines are plotted with formula (1); the black 
squares are the experimental results. The inaccuracy does 
not exceed 2%. A good agreement of theory and experi¬ 
ment is observed. The inset shows the residual difference 
between numerics (red lines) and calculations (black lines) 
with formula [TJ It is remarkable that at a relatively low 
magnetic field ~ 1000 G formula H]) describes the transi¬ 
tion frequency values with an inaccuracy of < 1 % (while at 
B = 1.8 kG the difference is < 0.3% which could be consid¬ 
ered as practically the full HPB regime). Fig. [TU] presents 


Fig. 8: Diagram of Di line transitions in HPB regime for 
laser excitation. The selection rules: Amj = 1; Ami = 0. 



Fig. 9: B-field dependence of the frequency shifts for transi¬ 
tions 1- 4 ot the Di line in case of a'^ excitation. Red lines 
are given by exact nnmerical calculations, black solid lines by 
formula ©, black squares by experiment. The inaccuracy is 
~ 2%. Inset presents the residual difference between numerics 
and calculations by formula (1) for curves 1-4- 



Fig. 10: The intensities and frequencies shifts (theory) of the 
transitions versus R-field for cr+ excitation. As seen for 
B > Bq only 1-4 transitions remain and the probabilities tend 
asymptotically to the same value at R 3> Bo- The inset shows 
F,mF for the ground and the upper levels for transitions 1-4 ■ 
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the theoretical values of the atomic transition probabili¬ 
ties (intensities) as a function of the S-field. As seen the 
intensities of transitions 1-4 at B > 1 kG tend to the same 
asymptotic value (HPB regime), which coincides with the 
experiment, while at B < 400 G, the intensity of the tran¬ 
sition labeled 4 (Fig. ITUl) is larger than those of 3 and 2 
which also coincides with the experiment. 



(a) (b) 


Fig. 11: a) Diagram of the hfs energy levels of Di line of in 
a magnetic field (for B < Bo) and possible Zeeman transitions 
for TT-polarized excitation are shown. The selection rules are 
AF — 0,±l;Am_F = 0. GT are labeled by rectangles; two 
IFFA transitions are indicated by dashed lines; b) Diagram of 
the K Di line transitions in the HPB regime for rr-polarised 
excitation. The selection rules are: Amj = 0; Am/ = 0. Two 
IFFA transitions are labeled by ovals. 

Magnetometry with A/2- and A- methods, tt- 
polarised excitation. — For this case the B-field is 
directed along the laser electric field E (magnets 3' are 
used and A/4 plate is removed). Fig. [Tl[ra)| shows the 
hfs energy levels diagram of the Di line of atoms in 
a magnetic field {B < Bq) and 14 possible atomic Zee- 
man transitions for rr-polarized excitation [including two 
forbidden at B = 0 atomic transitions mills] shown by 
dashed lines that we call ’’initially forbidden further al¬ 
lowed” (IFFA)]. Fig. inp^b)I presents the diagram of eight 
remaining transitions (the probabilities increase with the 
B-field). 

Note, that the necessary conditions for a modification 



Laser frequency detuning, MHz 


Fig. 12: Transmission spectrum of the NC with ®®K for L = A 
for B — 345, 420, 470 and 645 G ; vr-polarized excitation is used. 
GT and IFFA transitions are labeled by rectangles and ovals, 
respectively. The bottom curve is the reference. 



Fig. 13: B-field dependence of the frequency shifts for tran¬ 
sitions 1-8 of Di line for yr-polarized excitation. GT and 
IFFA transitions are marked by rectangles and ovals, respec¬ 
tively. The probabilities of six transitions (dashed lines) tend 
to zero for B > Bq. Solid lines - numerics; symbols - exper¬ 
iment (the inaccuracy is ~ 2%). The insets show the initial 
quantum numbers F, mp for the ground and upper levels. 

of the probability are the following nans]: the pertur¬ 
bation induced by the B-field couples only sublevels with 
mp — m'p = 0 which satisfy the selection rules AB = 0, 
AJ = 0, AB = 1, and modification of the probability is 
possible only for transitions between ground and excited 
sublevels when at least one of them is coupled with an¬ 
other transition’s sublevel according to the selection rules. 
As seen from Fig. ll]|la7| for the two side atomic transitions 
of K (marked as GT), |2,mF = +2) —>■ |2',mF' = +2) and 
|2,mF = —2) —?► |2',mF' = —2) the sublevels which could 
be mixed according to the selection rules are absent. Thus 
the probability of these two transitions remain the same 
in the whole range of applied B-fields, while the probabil¬ 
ities of other transitions differ significantly at low fields, 
but tend to the same value within the group at B ^ Bq. It 
is remarkable that the probabilities of all the other atomic 
transitions tend to that of two side atomic transitions (this 
is confirmed by the numerical calculations). For this rea¬ 
son why we call them guiding transitions. In Fig. [12] the 
transmission spectra of 7r-polarized excitation of the potas¬ 
sium NC (B = A for B = 345,420,470 and 645 G) are 
shown. GT and IFFA transitions are labeled by rectangles 
and ovals, respectively (as seen IFFA transitions undergo 
strong modification under the influence of the B-field). 

The B-field dependence of the frequency shifts for tran¬ 
sitions 1-8 for TT-polarized radiation is shown in Fig. 1131 
GT and IFFA transitions are labeled by rectangles and 
ovals, respectively. We see that transitions 1-8 are con¬ 
tained in two groups of four atomic transitions each. Note 
that the frequency slope (s) of four transitions inside the 
group asymptotically tends to the slope of GT. It is easy 
to show that they are equal to s = 0.94MHz/G and 
s = —0.94MHz/G for groups 1-4 and 5-8, respectively. 

The intensity and the frequency shifts (theory) of the K 
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Fig. 14: Theoretical curves for intensities and frequency shifts 
of the atomic transitions of Di line vs B-field for tt- 
polarized excitation. As seen for B ^ Bo only 1-8 transi¬ 
tions remain and the probabilities tend asymptotically to the 
same value as for GT transitions (by module 4.66x 10~^® ESU). 
Groups 1-4 and 5-8 are located at the high and low frequency 
side, respectively, while the transitions labeled 1 and 8 have the 
largest and smallest frequencies, respectively. GT are marked 
by arrows (each by the same colour) and as seen their ampli¬ 
tudes remain the same in the range of 100 — 800 G. 

transitions vs B-field for 7r-polarized excitation are shown 
in Fig.[T4l As seen for B ^ Bq only transitions 1-8 remain 
and the probabilities at i? 3> So tend asymptotically to 
the same value as GT (by module 4.66 x 10“^® ESU). The 
inset shows the diagram of the K Di line transitions in 
HPB regime with GT and IFFA marked. 

Summary and conclusions. — It is remarkable that 
with K atoms (due to the small value of Bq for ®®K) con¬ 
tained in NG all the peculiarities of the behavior of atomic 
transitions of the Gs, ®’’Rb, and ®®Rb in the whole region 
of S-field up to the HPB regime could be easily studied 
using smaller, at least by an order, S-fields. For example, 
in Ref. [53] the ®^Rb Di line for TT-polarized excitation 
was studied (the diagram is similar to the one shown in 
Fig. ima)). For a R-field of 220 G there was no any ev¬ 
idence of IFFA transitions, while for the K atoms IFFA 
transitions (shown in Fig. 1121) have amplitudes compara¬ 
ble with the other transitions. Also, the peculiarities of the 
HPB regime for the K atoms can be detected at H-fields 
smaller by at least one order. Note that in [T] a magne¬ 
tometer based on the Rb NG use was described, while the 
spatial resolution of K NC for L = A/2 = 385 nm should 
be better than that of the Rb (398 nm) and Cs (448 nm), 
which is important when measuring a strongly inhomo¬ 
geneous magnetic field. It is worth to note that in some 
cases instead of Rb NG a micrometric-thin (^ 40 ^m) cell 
could be successfully used m- Thus, we suppose that 


with the help of a 40 p,m-thin K cell the above mentioned 
peculiarities could be studied, too. 

* * * 

The authors thank A. Papoyan, G. Adams, L Hughes 
and H. R. Jauslin for discussions. The research was con¬ 
ducted in the scope of the International Associated Labo¬ 
ratory IRMAS (GNRS-France & SGS-Armenia). A.S. and 
D.S note that this work was also supported by State Gom- 
mittee Science MESRA, in frame of the research project 
No. SGS 13-1C029. A.S. and A.T. thank the ANSEF 
Fund - Opt 3700 grant. 

REFERENCES 

[1] Sargsyan A., Hakhumyan G., Papoyan A. et al, Appl. 
Phys. Lett, 93 (021119) 2008. 

[2] Tremblay P. et al, Phys. Rev. A, 42 (2766) 1990. 

[3] Budker D., Gawlik W., Kimball D. et al. Rev. Mod. 
Phys., 74 (1153) 2002 

[4] Auzinsh M. et el. Optically Polarized Atoms: Under¬ 
standing Light-Atom Interactions, edited by , Vol. (Oxford 
University Press, Oxford, New York) 2010, p. 55-62. 

[5] Bloch D., Ducloy M. et al. Laser Physics, 6 (670) 1996 

[6] Umfer C., Windholz L., Musso M., Z. Phys. D: At. 
Mol. Clusters, 25 (23) 1992 

[7] Olsen B.A. et al, Phys. Rev. A, 84 (063410) 2011 

[8] Sargsyan A., Hakhumyan G., Leroy C. et al. Opt 
Lett , 37 (1379) 2012 

[9] Weller L., Kleinbach K.S., Zentile M.A. et al. Opt. 
Lett, 37 (3405) 2012 

[10] Sargsyan A. et el, JETP Lett, 98 (441) 2013 

[11] Sarkisyan D., Bloch D., Papoyan A., Ducloy M., 
Opt. Commun., 200 (201) 2001 

[12] Keaveney J., Sargsyan A., Krohn U. et al, Phys.Rev. 
Lett, 108 (173601) 2012 

[13] Briaudeau S., Saltiel S., Nienhuis G., Bloch D., 
AND Ducloy M., Phys. Rev. A, 57 (R3169) 1998 

[14] Dutier G. Yarovitski A., Saltiel S. et al, EPL, 63 
(35) 2003 

[15] Sarkisyan D., Varzhapetyan T., Sarkisyan A. et al, 
Phys. Rev. A, 69 (065802) 2004 

[16] Pashayan-Leroy Y., Leroy C., Sargsyan A. et al, 
J. Opt Soc. Am. B, 24 (1828) 2007 

[17] Zentile M.A., Andrews R., Weller L. et al, J. Phys. 
B: At Mol. Opt Phys, 47 (075005) 2014 

[18] Hakhumyan G., Leroy C., Pashayan-Leroy Y. et al. 
Opt Commun., 284 (4007 ) 2011 

[19] Sargsyan A., Hakhumyan G., Leroy C. et al, ,J0SA 
B, 31 (1046) 2014 

[20] Hakhumyan G., Leroy C., Mirzoyan R. et al, EPJD, 
66 (119) 2012 

[21] Zentile M. et al. Computer Physics Commun., 189 
(162) 2015 

[22] Sargsyan A., Tonoyan A., Hakhumyan G. et al. 
Laser Phys. Lett, 11 (055701) 2014 

[23] Sarkisyan D. et al, J. Opt. Soc. Am. B, 22 (88) 2005 

[24] Sargsyan A. et al. Opt. Lett., 39 (2270) 2014 


p-6 



































